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Abstract. We use iV-body simulations to find the effect of neutrino masses on halo 
properties, and investigate how the density profiles of both the neutrino and the 
dark matter components change as a function of the neutrino mass. We compare 
our neutrino density profiles with results from the A^-one-body method and find good 
agreement. We also show and explain why the Tremaine-Gunn bound for the neutrinos 
is not saturated. Finally we study how the halo mass function changes as a function 
of the neutrino mass and compare our results with the Sheth-Tormen semi-analytic 
formulae. Our results are important for surveys which aim at probing cosmological 
parameters using clusters, as well as future experiments aiming at measuring the cosmic 
neutrino background directly. 



1. Introduction 

Massive neutrinos are known to have a significant effect on cosmic structure formation [H 
[2] . In the early universe they contribute to the relativistic energy density and influence 
the transition from radiation to matter domination. At late times they contribute to 
the dark matter density, and therefore also to cosmic structure formation. However, as 
opposed to Cold Dark Matter (CDM), they do not contribute to structure formation 
on physical scales smaller than the free-streaming scale, roughly equal to the distance 
traversed before the neutrinos become non-relativistic. This suppression of small-scale 
structure leaves a very distinct imprint on large-scale structure observables such as the 
matter power spectrum, which can in turn be used to probe neutrino physics. Many 
studies have by now been devoted to this topic, most of which focussing on constraining 
the neutrino mass, m^. At present an upper limit on the neutrino mass can be derived 
from observations of the Cosmic Microwave Background (CMB) anisotropics alone or 
in conjunction with various large-scale structure data sets, such as the SDSS-DR7 LRG 
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catalog, and falls in the range ^ m^, < 0.4 — 0.7 eV, depending both on the complexity 
of the model space and the combination of data sets used (e.g., [21 IH IS])- In the future 
the sensitivity of large-scale structure observations to the neutrino mass will increase 
significantly. For example it has been estimated that the combination of CMB data 
from Planck and a weak lensing survey from the LSST will push the la sensitivity to 
better than 0.05 eV, close to the minimum Yl'^f allowed by oscillation data [6]. 

While most neutrino mass constraints at present have been derived using large-scale 
structure correlation functions (or power spectra), there are other observables that are 
potentially just as interesting. One prime example is cluster number counts which are 
in principle very sensitive to the neutrino mass [7]. However, in order to fully utilise 
such data it is necessary to have accurate theoretical predictions, which so far do not 
exist for ACDM models extended with massive neutrinos (see, however, [8] for an early 
calculation based on the old mixed dark matter scenario). In the present paper we 
calculate the halo mass function in ACDM cosmologies with massive neutrinos included 
for a variety of neutrino masses. However, before proceeding to this and a discussion of 
other observables related to halo properties, let us briefly review how neutrinos affect 
structure formation in the linear regime. 

1.1. The effect of neutrinos 

The effect of neutrinos on structure formation in linear theory has been studied numerous 
times in the literature (see, e.g., [9]). In general the amount of fluctuations at a given 
wavenumber k is represented by the power spectrum, P{k) = which can be split 
in the following form 



where D{z) is a scale independent growth factor and Pq is the initial power spectrum. 
T{k) is the transfer function (TF) which is both time and scale dependent in general. 
The effect of massive neutrinos is embedded entirely in the TF and is separated into 
two regimes. On scales much larger than the free-streaming scale. 



where rui, is the one-particle neutrino mass, neutrinos behave essentially like CDM, while 
on smaller scales they suppress structure formation. Very naively one might expect the 
suppression arising from replacing a fraction of the CDM component with neutrinos to 
be of order AP/P ~ —Q^/Qm, where Qrn = ^Ic + + is the total matter density, 
because neutrinos do not cluster. However, this grossly underestimates the true effect 
because massive neutrinos also influence the background expansion around the time of 
matter-radiation equality. The final result in linear theory is that the suppression is 
approximately given by 




(1) 



A;fs ~ 0.8— y/iMpc 



eV 
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This shows that most of the effect actually comes from the modification to the 
background, i.e., sub-eV to eV scale neutrinos lead to a longer radiation era. This 
effect is also much larger than the effect of replacing a fraction AQrn of the CDM energy 
density with A. In this case the change in the matter power spectrum on small scales 
(with the large-scale normalisation held constant, i.e., ignoring the effects of AQm on 
the growth factor) is very approximately given by 



approximately two times smaller than that due to assigning AQm to massive neutrinos. 

Since neutrinos have such a strong effect on the power spectrum even in linear 
theory it is natural to expect a similarly strong effect in the non-linear regime. This 
was tested in detail for the power spectrum in a number of papers (TDl [H] [T2], and 
a significant enhancement in the power spectrum suppression was indeed found: The 
maximum suppression is increased from to approximately —9.8^2,^/17^13 with 

a pronounced feature at ~ 0.7/iMpc"^. 

Another issue which has so far not been addressed with precision A^-body 
simulations is how the presence of massive neutrinos affect halo formation. Here we 
study how CDM halo properties are altered by the presence of massive neutrinos, and 
we also present detailed results for the corresponding neutrino halos. This last point 
is important for example for understanding the prospects for a direct experimental 
detection of the cosmic relic neutrino background. 

The paper is organised as follows: In Section [2] we present the numerical setup 
required for the analysis. In Section [3] we present results on halo profiles for both the 
neutrino and matter components. In Section H] we discuss how the halo mass function is 
altered in models with massive neutrinos, and finally Section [5] contains our conclusions. 

2. Numerical setup 

2.1. Initial conditions and N-body simulations 

Our A^-body simulations are carried out using the modified version of the gadget-2 
code [H] described in [12]. The code utilises a hybrid scheme for simulating the neutrino 
component: All neutrinos with velocities much higher than the average gravitational 
flow velocities in the simulation are treated using linear perturbation theory, and their 
effect on the gravitational potential in the A^-body simulation is included via a Fourier 
grid. At low redshift the low velocity neutrinos, with q/T < 6, are followed separately 
in 6 momentum bins by A^-body particles {q is the comoving momentum and T the 
comoving temperature). In Fig. |4] it can be seen that the higher part of momentum 
space, q/T > 6, do not contribute to neutrino clustering in halos. 

I This finding was confirmed in tlie very recent paper [13| . 




(4) 



where 0,1 



'm 



Qm + Afim, for small changes in AQrn ^ Qm- The effect here is 
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The initial conditions (ICs) for the CDM and neutrino components are generated 
with the same set of random numbers for a given box size. This reflects the assumption 
of adiabatic primordial ICs. The CDM distribution is followed with CAMB [15] until 
z = 49, where the Zel'dovich approximation [TB] and a second-order correction calculated 
with second-order Lagrangian perturbation theory [17] are used to generate the CDM 
iV-body particle initial displacements and gravitational flow velocities. 

When the maximum thermal velocity of the first neutrino momentum bin, q/T = 1 
(upper limit), has fallen below /flow = 4 times the average CDM gravitational velocityj§| 
the low momentum neutrinos, which have non-linear clustering, are converted to A^- 
body particles with the Zel'dovich approximation and followed with gadget-2. The 
high velocity part is retained on the grid. The ICs for each bin are generated 
from a momentum dependent TF, and in addition to the gravitational flow velocities 
the neutrino A^-body particles receive a thermal, Fermi-Dirac distributed, velocity 
corresponding to the particular bin. In effect we are simulating 6 different Hot Dark 
Matter components with different thermal properties. We employ a new timestep 
criterion for the neutrino component oc {1 + z)/Q^. 

The TF for each momentum bin is calculated with CAMB. By default CAMB 
uses 15 neutrino momentum bins, which is sufficient for sub-percent accuracy in the 
momentum averaged TFs. However, for our purpose, we need percent level accuracy 
in each individual momentum bin, and to ensure this we use a total of 480 individual 
TFs such that each of the 15 momentum bins is reconstructed from 32 separate TFs. 
See [12] for further information on the hybrid implementation of neutrinos in A^-body 
simulations. 

The hybrid method has the advantage that all neutrino masses can be simulated 
without compromising computational speed or accuracy (treating low mass neutrinos 
as particles leads to prohibitively high CPU time consumptioiJJJ] and particle shot noise, 
and, conversely, treating high mass neutrinos in linear perturbation theory leads to 
loss of accuracy). The code allows us to accurately calculate halo number densities 
and properties over a wide mass range. Since neutrinos cluster much less than CDM 
ultrahigh resolution is in most cases not necessary]^ For example, since the neutrino halo 
profiles flatten at small r one does not encounter the same resolution problem as with 
CDM. Likewise there is hardly any neutrino halo substructure because neutrino halos 

§ When the box size is increased the average CDM gravitational velocity increases as well due to 
extra large-scale velocity flows. As a result, the neutrino A^-body particles are created earlier in 
large simulation volumes, and not due to larger non-linearities which of course would justify a larger 
conversion rcdshift. But this effect is small, at the order of 5% in redshift. Ideally, the velocities should 
be interpolated to a grid, Fourier transformed and then convolved with a window function eliminating 
large-scale velocity modes. 

II Further investigations for an appropriate Courant timestep condition for the neutrino component 
could alleviate this problem. 

% Ultrahigh resolution can be achieved by making zoom simulations, though for neutrinos the advantage 
of such a setup is smaller, since fast-moving neutrinos can easily leak out of the high-resolution region, 
and hence the initial volume of the zoomed neutrinos has to be considerably larger than the initial CDM 
volume (though a periodic volume around the halos could alleviate this problem at higher redshift). 
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Table 1. A^-body simulation parameters: iVcDM and A^,y,part are the number of CDM 
and neutrino A^-body particles respectively, A'^^.grid the size of the linear neutrino 
Fourier grid, Qcut/T the cut-off below which the neutrino component is converted to 
particles, and /flow determines the redshift of this conversion. The simulation box 
size is represented by i?BOX, ''^f the total neutrino mass roughly related to the 
neutrino density parameter, fij^, hy ft^ = m^/ (94 h^cY), with Qm = + fib + fl^,. 
The exotic simulations E1-3 have no 61, in the iV-body simulation, but the neutrinos 
are still included in the background evolution, while £'4-5 have a CDM iV-body particle 
mass corresponding to 0,„ = 0.3 and the ICs are calculated from weighed CDM and 
baryon TFs from a cosmology with neutrinos included. 

are primarily formed by late time neutrino infall on already existing CDM structures. 

Parameter setup We assume a standard fiat cosmology with h = 0.7, as = 0.878 (for a 
model without massive neutrinos), = 1, = 2.3 ■ 10~^(r2m/0.3)^, Qf, = 0.05 and 
varying amounts of Qc and In cosmologies with massive neutrinos we assume 

3 degenerate neutrino species. Table [1] shows parameters for the various A^-body 
simulations presented in this paper. 

By running simulations in simulation volumes of 32, 128 and 256 Mpc we found 
that the neutrino halo profiles were almost identical in the latter two box sizes, and that 
the smallest box significantly affected the density profiles for the larger halos. In sum, 
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we only present results for simulation volumes larger than 256 Mpc. 

We have chosen /aqw = 4 as the criterion for creating neutrino A^-body particles, 
which we have shown to be a reasonable value by comparing with simulations in which 
/flow = 2 and 8. 

Finally, when we present matter density profiles for cosmologies with = 
0.15 eV and 0.3 eV neutrinos, only the homogeneous neutrino component has been 
added. This seems reasonable since these low mass particles contribute insignificantly 
to the overall matter density profile. 

Fig.[T]shows the CDM and neutrino A^-body particles in halos with different masses. 
Here, individual A^-body particles can be identified, with small and bright particles lying 
in high density regions, and larger and darker particles in lower density areas. 

Fig. [1] clearly illustrates the lack of neutrino particle statistics in halo centers. 
Neutrino clustering is determined by the combined effects of gravity and thermal 
velocity. The smaller the neutrino and the halo masses the less clustering. To simulate 
such small overdensities a very fine initial A^-body particle grid is required. From the 
figure it can also be seen that only q/T < 3 trace the underlying CDM distribution on 
the scales shown, whereas structure in the higher momentum bins can only be seen on 
larger scales. This is also consistent with Fig. HI 

2.2. Halo finding 

We use the Amiga Halo Finder (AHF) [18] to identify halos and their centers within the 
simulations. Only bound particles have been used to calculate halo centers and their 
virial radii. From the halo centers we calculate the matter and the neutrino halo profiles 
for halo masses of 10^^, 10^^, 10^^ and 10^^ Mq. We show matter halos and not CDM 
halos, since it is the former quantity which is the measurable one. 

We stack halos in mass bins with a bin width of 10% of the central halo mass. These 
widths are narrow enough to ensure a reliable halo profile calculation of the desired halo 
masses, while at the same time providing enough halos to stack. The A^-body density 
profiles presented in this paper have been found by fitting a smooth curve through the 
data points, though we note that this is mainly important in the inner part of the 
neutrino density profiles for low neutrino and halo masses, where particle statistics is 
low. 

Furthermore, to compare with the A^-one-body method, we have only used halos 
with an average overdensity within the virial radius, Ayir, in the range 330 — 340 for 
Qm = 0.3 (340 — 350 for Q^n = 0.287). This criterion eliminates, e.g., gravitationally 
stripped halos with a very high central density. 

2.3. The N -one-body method 

The A^-one-body method was introduced in [19]. It is a restricted method devised to 
solve, approximately, the (non-linear) collisionless Boltzmann equation for the neutrino 
phase space based on the following observation: In the limit pu ^ pm density 
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Figure 1. CDM and ^mi, = 1.2 eV neutrino distributions for halo masses ~ 
5 • IO^-^Mq (top), ~ IO^-^Mq (middle) and ~ IQ^'^Mq (bottom), where the masses 
only correspond to the central halos in the upper two mosaics. Dimensions in each 
image are 5, 2 and 1 Mpc, respectively. In each mosaic the images correspond 
to CDM, total neutrino, and q/T = 1 to 6 from top-left to bottom-right. Individual 
neutrino A'^-body particles can be identified. 
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perturbations in the CDM fluid dominate the total gravitational potential, and not only 
will the CDM halo be gravitationally blind to the neutrinos, gravitational interaction 
between the neutrinos themselves will also be negligible. This allows us to treat the CDM 
halo as an external gravitational source, and compute the trajectory of each neutrino 
phase space element as a test particle moving in an external potential one at a time in 
N independent simulations. An obvious advantage of this technique is that it requires 
virtually no computing power when compared with a full scale A^-body simulation with 
the same, large A^, and is thus particularly useful for resolving the clustering of neutrinos 
on small scales. 

In principle the A^-one-body method can be applied to any given CDM density 
distribution. Here, however, we focus exclusively on spherically symmetric CDM halos 
whose density profiles are parameterised by fitting functions. Apart from possible 
deviations in the innermost cores, dark matter halos are well described by the universal 
Navarro- Frenk- White (NFW) profile [20], with a density given by 

(r/rs)(l + r/rs)^ 

where = 4p(rs) is a characteristic inner density at the characteristic radius rg. r^ is 
related to the virial radius rvir via the concentration parameter c 

c=^-^. (6) 

In a ACDM cosmology, a good analytical expression for the concentration parameter 
is [21] 

where M^^^ is the virial mass lying within rvir, and c{z) ~ c(0) /{l + z). We expect c to be 

smaller for cosmologies with neutrinos since the neutrinos free-stream out of the inner 

density cores and effect the halo formation process. Defining the average overdensity 

within Tvir to be A^ii., we have 

4?? r c 1 

Mvir = — AvirPm.or-vir = 4:npsa^r^ ln(l + c) - — — , (8) 

O LI ~\- C- 

so that a halo of a given mass is fully described in terms of its concentration parameter 
alone. The overdensity Avir can also be approximated by the overdensity at virialisation 
from the spherical top-hat collapse model, 

. 187r^ + 82[nmiz) - 1] - 39[n^{z) - 1]^ 

with ^lm{z) = fim/(^m. + ^^A^^) [22] • For = 0.3, we find 5th — 337 (5th — 346 for 
Qrn = 0.287), which explains why we only use the range Avir = 330 — 340 (340 — 350) 
when analysing A^-body data. 

We model the CDM distribution as a NFW halo sitting on top of a uniform 
distribution of CDM, i.e., the A^-one-body method assumes that all halos are completely 
isolated. In order that the halo overdensity merges smoothly into the background 
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density, we extend the NFW profile to beyond the virial radius. The initial neutrino 
distribution is taken to be spatially uniform0 with a momentum distribution described 
by relativistic Fermi-Dirac statistics. This initial distribution is divided into small 
chunks in both real and momentum space, and each chunk is allowed to move under the 
external potential of the CDM halo, but independently of each other. A low resolution 
run is first carried out for each set of neutrino and halo masses. All chunks that end 
up at z = inside a sphere of radius 50 Mpc centered on the halo are traced back 
to their origin, subdivided into smaller chunks, and then re-simulated. The process is 
repeated until the inner ~ 10 kpc is resolved. 

The A^-one-body simulations are started at the same redshifts as when the neutrino 
A^-body particles are created. This redshift is neutrino mass and halo mass dependent. 
Note that when comparing the density profiles with the two methods, the underlying 
CDM profiles are not identical, since the neutrino component do contribute to the 
gravitational potential and not least the halo merger history in the A^-body simulations. 

3. Halo structure 

3.1. Neutrino clustering 

The gravitational effect of a host halo is relatively much more important for neutrinos 
than for the CDM component: Due to free-streaming neutrinos will almost completely 
stream out of small halos (~ IO^^Mq), and any measured value 6^ > will be caused by 
the host halo. The radial profile of 6^, will therefore be a superposition of a dominant 
flat profile from the host halo on top of a sub-dominant contribution from the ~ lO^^M© 
halo itself. This fact can be seen in Fig. [21 

Since it is the isolated halo profile that the A^-one-body method calculates, we 
have found the neutrino density profiles from isolated halos in the A^-body simulations. 
We used the criterion that for a halo to be considered isolated it should be more 
than 10 times the virial radius away from a heavier halo. With this criterion the 
agreement between the two methods improves significantly for masses < lO^'^M©, 
in effect confirming the robustness of both the A^-body and the A^-one-body methods. 
It can also be seen that for = 0.3 eV there is a significant contribution to 6^, even 

in a ^ IO^^Mq halo from heavier halos, which is not the case for more massive neutrino 
states. This is caused by the fact that as the neutrino mass decreases, neutrinos free- 
stream out of ever larger halos, so that the relative effect of even larger host halos must 
be taken into account. 

The effect of tidal truncation on the A^-body halo profiles can easily be seen in 
Fig. [2j This effect is not included in the A^-one-body approach. Furthermore, for the 
IO^^Mq halos it can be seen that the neutrino density falls below its cosmic average 
beyond the virial radius. This could be due to either the presence of underdensities at 

We have run one A'^-body simulation where the neutrino particles were assigned homogeneous ICs, 
which in general confirmed the validity of this assumption. 
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Figure 2. Neutrino halo profiles for ^rrii, — 0.3 eV (top), ^rrii, ~ 0.6 eV (middle) 
and J2fn^ = 1.2 eV (bottom) for halo masses of 10^^ lO", lO" and 10^^ Mq. Profiles 
are calculated with the iV-one-body method (dotted) and the A^-body method with a 
halo isolation criterion (solid) and without (dot-dashed). 
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Figure 3. Halo profiles from iV-body simulations for a model without massive 
neutrinos, with isolated halos (solid) and all halos (dot-dashed). The halo masses 
are 10^^, 10^"^, 10"'^"' and 10^^ M©. The profiles for the lowest 3 halo masses are taken 
from the 256 Mpc box and the profile for the most massive halo is taken from the 
1024 Mpc box. NFW profiles are also shown (dotted), and the halo mass dependent 
virial radii are indicated by the '+' signs. 



particular distances from the halo centers, or due to the fact that we only select isolated 
halos, which are more likely to be found in low density regions. 

From the pure ACDM A^-body simulations presented in Fig. |3] it can be seen 
that our matter halos are perfectly fitted by a NFW profile over the mass range 
Mvir = 10^^ — IO^^Mq until 20 kpc from the halo centers. Here our A^-body results 
begin to lack particle resolution. The profile for the larger halo mass is taken from a 
1024 h-^ Mpc box with the same number of particles, and this halo is therefore only 
resolved until ~ 100 /i~^kpc. Note that our dominant background NFW profiles in 
the A^-body simulation are valid down to scales significantly smaller than the scales at 
which we present neutrino density profiles. Therefore, our neutrino density profiles are 
not affected by insufficient CDM A^-body particle resolution. 

Since the CDM component is much more clustered than its neutrino counterpart, 
the flat profile from the host halo is only dominant relative to the contribution from 
the halo itself on scales beyond the virial radius (see Fig. From this figure it can 
also be readily understood why the neutrino density profiles differ when only low mass 
isolated halos are considered: The underlying CDM gravitational source term is roughly 
fiat beyond the virial radius, and within the virial radius the neutrinos free-stream out 
of the small mass halos, in sum producing a roughly fiat neutrino density profile also 
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Figure 4. Cumulative neutrino halo density profiles as a function of momentum in 
a IO^^Mq halo, for J^m^ = 0.6 cV (left) and Y^m^ = 1.2 cV (right). The neutrino 
density with q/T > 6 has been added as a homogeneous term to all profiles. 



within the virial radius. 

In Fig. Hlwe show the cumulative neutrino density profile for a lO^^M© halo, for total 
neutrino masses of 0.6 eV (left) and 1.2 eV (right). It can be seen that only neutrinos 
with q/T < 3 contribute within the inner 1000 kpc, whereas neutrinos with momenta 
up to q/T ^ 5 is needed to simulate profiles beyond ~ 2000 kpc. It can also be seen 
that as the neutrino mass is increased higher neutrino momentum bins contribute to 
neutrino clustering on a given scale. Finally, the figures confirm the accuracy of only 
converting neutrinos with g/T < 6 to A^-body particles. Since this result is accurate for 
a 1O^^M0 halo, it is certainly also accurate for lower mass halos where only neutrinos 
with a momentum from the very low end of the Fermi-Dirac distribution cluster. 

The Tremaine-Gunn bound Based on purely theoretical grounds one should expect the 
following scenario. Neutrinos cluster in a halo of given mass M^i,. and radius r^^,., so 
that the escape velocity at rvir of a given halo is 



.^y^. (10, 

V ' vir 

Only neutrinos up to this velocity can be bound in the halo and the central density 
should therefore be 

myp^f{p)dp. (11) 
For small p, / is approximately 1/2, and we have 

P^-{^T<. (12) 

\ 'vir / 



leading to 



6,^[^\"\l. (13) 



vir 
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Figure 5. The evolution of the neutrino density profile around a static NFW halo of 
IO^'^Mq for '}2ra^ = 0.3 eV (left) and "Y^mi, = 0.6 eV (right). The physical density 
perturbation, a~^5^, as a function of the physical radius, ar, is time- independent in 
the TG limit. 



This is the Tremaine-Gunn (TG) bound [23], and is essentially the same result quoted 
in Eq. (7.3) of [l9] for a NFW halo. 

However, in practice the neutrino halo density almost never saturates the Tremaine- 
Gunn bound, as was also seen in Fig. 7 of [19] . Fig. |5] sheds some light on the actual 
evolution of the neutrino halos. The setting is a simplified version of the A^-one-body 
method in which the NFW halo is taken to be static so that the TG bound is time- 
independent in physical units. For the smallest neutrino mass ^ mj, = 0.3 eV, it can 
be seen that the physical density perturbations of neutrinos drop as the scale factor 
increases, stabilising only at a ~ 4. At early times, i.e., a = 1/2, the TG bound is 
almost saturated because the neutrino density contrast is low and the halo is populated 
only with neutrinos drawn from the very low momentum end of the relativistic Fermi- 
Dirac distribution for which / ~ 1/2. However, this also means that the final neutrino 
phase space density deviates little from the relativistic Fermi-Dirac distribution, so that 
linear perturbation theory remains valid. The dilution due to the background expansion 
dominates until very late and the physical neutrino density perturbation only approaches 
a constant at a ~ 4, and at a value much lower than the TG bound. 

For the higher neutrino mass this effect is less pronounced, i.e., the physical density 
levels off much earlier. However, for the higher neutrino mass the TG bound is far from 
saturated even at early times. The reason for this can be understood from Fig. HI For 
the higher neutrino mass, neutrinos with q/T > 1 make up most of the halo. However, 
since the initial neutrino phase space density is far less than 1/2 at q/T > 1, it also 
makes it more difficult for the final coarse-grained phase space density to reach the 
upper limit of / ~ 1/2 at the higher momentum end of the spectrum. For example, the 
second bin with 1 < q/T < 2 has / = P{q)q'^dq/ f\q)q^dq ~ 0.18. 
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Figure 6. Change in the matter halo profiles relative to a base model without massive 
neutrinos for a halo mass of 10^^ M0. Left: As a function of neutrino mass. Right: 
For different (exotic) cosmologies (see details in text). 



3.2. Feedback on CDM halos 

Even though neutrinos make up a minute fraction of the dark matter halo mass at late 
times, they do have an impact on the halo formation history, i.e., they affect the halo 
mass dependent merger rate [21]. In this section we are interested in alterations of the 
halo profiles, not the number density of halos in a given mass interval (i.e., the halo 
mass function). This aspect will be discussed in detail in Section |H 

Fig- El (left) shows the relative change of the matter halo density profile as a function 
of neutrino mass for a 10^^ M0 halo. We only show results for this halo mass, since here 
the product of the number density of halos times the halo density is maximal. From 
the figure it can be seen that the presence of massive neutrinos lowers the density at 
radii smaller than ~ 100 kpc. Since we compare halos of identical total masses this 
is compensated by an increased density at larger radii. The reason for this effect is 
that halos form later in models with massive neutrinos because the linear theory TF is 
lowered. From pure CDM simulations it is indeed known that late forming halos are 
less concentrated [25], with the concentration parameter scaling roughly as c oc l/oc, 
where is the scale factor at formation. 

From the right hand side of Fig. [6] it can be seen that removing the neutrino 
perturbations (labelled 6y = 0) from a cosmology with 'Y^rriu = 0.6 eV in the A^-body 
simulation leads to less concentrated halos. On the scales shown the suppresion is as 
large as 3 — 4%, though we caution that there is some noise in the data. It can also 
be seen that replacing Q^, by an enlarged Qc (labelled Qc = 0.3) when the A^-body 
simulation is started gives a profile which is too steep at the few % level. Finally, a 
cosmology with Qrn = 0.287 (labelled = 0.287) is placed somewhere between the 
aforementioned cosmologies. 
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4. The halo mass function 

We now turn our attention to the number density of halos in a given mass interval, the 
halo mass function (HMF). In Fig. [7] (top and middle figures) we show the HMFs for 
cosmologies with different neutrino masses. As expected the HMF is more suppressed 
in cosmologies with a larger neutrino mass. The suppression is largest for the heaviest, 
late forming halos. 

The last panel in Fig. [7] shows that the cosmologies where the Yl'^^f — 0.6 eV 
neutrino component is replaced by CDM at the redshift where the A^-body simulation 
is started as well as the one with Qm = 0.287 predict HMFs with different shapes. Only 
a model without neutrino perturbations in the A^-body simulation gives roughly the 
same result as the full calculation: The suppression of the HMF is mainly caused by the 
suppression of the initial TF in the linear regime and not by neutrino clustering effects 
in the A^-body simulation. 

In Fig. [8] we compare our HMFs calculated from A^-body simulations with HMFs 
from the Sheth-Tormen (ST) semi-analytic formulae [2S]0 The ST fit is based on the 
fact that, as first pointed out by Press and Schechter [29], the HMF can be written as 

MdM dn(M,z) ,du , 

—p dM- = '^^''^V^ 

with u = [Ssciz) / a{M)]'^ , where 6sc{z) = 1.686 is the overdensity required for spherical 
collapse at z, and p = ^ImPc- dn{M, z) is the number density of halos in the mass 
interval M to M + dM. The variance of the linear theory density field, o"^(M), is given 
by 

<y\M) = J f^^^mkR)\\ (15) 

where P\in{k) is the linear theory matter power spectrum, and the Top-Hat window 
function is given by W{x) = (3/x^)(sina; — xcosx) with R = {SM/AnpY^^. 
The ST fit to uf{iy) is 

+ (16) 



with u' = 0.707z/ and p = 0.3. A = 0.3222 is determined from the integral constraint 
Jf{u)du=l. 

The upper panel in Fig. |8] shows that the agreement is poor if flm = ^c + ^b + ^u is 
used in the ST formalism. However, this is due to a wrong definition of the halo mass: 
Even for the very largest cluster halos the neutrino component contributes very little to 
the halo mass. In reality, the mass inside the collapsing region should be calculated using 
Qc + ^b, not Qm- This amounts to neglecting the weakly clustering neutrino component 
when calculating the halo mass. The two lower panels in Fig. [S] shows the same ST fit, 
but using Qc + ^6 instead of Qm- In this case the ST HMFs provide an excellent fit 

* See [17J for other semi-analytic fits to the mass function, which we found to work less well for 
models with massive neutrinos. 
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Figure 8. Relative halo mass functions for different neutrino cosmologies compared 
with the predictions from the Sheth-Tormen formulae (black lines). Top: With 
ilm = + + m the ST formulae. Middle: With Oc + f^b used instead of in 
the ST formulae. Bottom: Differences between the iV-body and the ST predictions. 
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to the relative change to the HMF caused by neutrinos. As the figure at the bottom 
clearly demonstrates, the agreement is better than ~ 3% at halo mass scales where our 
A^-body HMFs are accurate (we do not consider our A^-body HMFs to be accurate by 
more than a few % on any mass scale, M). Although the absolute HMFs, even for CDM 
simulations, do not match the ST HMFs more precisely than at the ~ 10% level, the 
relative change from adding neutrinos can be calculated significantly more accurately. 

5. Conclusion 

We have performed a detailed study of halo properties in ACDM cosmologies with 
massive neutrinos included. An important goal was to study the neutrino density 
profiles in dark matter halos. To this end we employed detailed A^-body simulations 
across a wide range of scales to test halo masses from Milky Way size (10^^ Mq) to 
large clusters (10^^ M©), as well as the A^-one-body method developed to solve the 
neutrino Boltzmann equation approximately around existing CDM halos. In general 
we found good agreement between the full A^-body and the A^-one-body results. The 
difference between the A^-body and A^-one-body methods arise from the fact that the 
latter assumes the CDM halo to be monolithic and at all times describable in terms of 
a NFW profile, i.e., it does not take into account halo substructure and larger merger 
events. It also assumes an analytic evolution of the concentration parameter. We also 
discussed in some detail how the density profiles of neutrino halos can be understood in 
terms of the Tremaine-Gunn bound, i.e., the bound coming from the fact that a coarse- 
grained distribution can never attain values exceeding the maximum of the original 
fine-grained distribution. 

For smaller halo masses, the neutrino profiles in isolated halos are in excellent 
agreement with the prediction from the A^-one-body method. This result is not too 
surprising since this is exactly the case where the infall on an existing spherical NFW 
halo is most realistic. However, many smaller mass halos are embedded in larger cluster 
halos and for the smaller neutrino masses the local neutrino profile in such a halo is 
dominated by the background of neutrinos bound in the much larger cluster halos. 

In terms of the local neutrino density enhancement, which is relevant for possible 
future attempts at direct Cz/B detection, a Milky Way-size galaxy halo is too small 
to have a significant overdensity, even when taking a possible cluster background into 
account. 

We also briefly studied how neutrinos impact on the density profiles of the CDM 
halos. While neutrinos contribute very little to the total density in the halo, the presence 
of massive neutrinos in the model leads to slightly later formation of halos with a given 
mass and consequently to generally lower concentration parameters, c. 

Finally, we calculated halo mass functions for ACDM models with massive 
neutrinos. Since large cluster surveys will become available in the coming years, the 
halo mass function is an important cosmological observable. As expected, we find a 
very strong suppression of halo formation with increasing neutrino mass. As noted in 
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previous analytic or semi- analytic studies the suppression is particularly marked for 
massive halos because the suppression in linear theory power from massive neutrinos 
shifts the maximum cluster mass down, i.e. the scale beyond which the halo mass 
function is exponentially suppressed. 

We then compared the halo mass functions from simulations with halo mass 
functions calculated using the semi-analytic method developed by Sheth and Tormen. 
If used naively, i.e. just processing the linear theory power spectrum without any 
adjustment to the method, the agreement is poor. However, it is easy to see that 
the disagreement arises because the ST method implicitly assumes that all matter 
clusters in the same way (the value of Qm used is Qc + + ^u)- However, even 
large clusters bind relatively few neutrinos and for all halos it is true that neutrinos 
make a negligible contribution to the halo mass. If the ST formalism is corrected for 
this by using Qrn = + ^b, i-e. taking into account only the clustering species (but of 
course using the correct initial power spectrum and the correct background evolution), 
the agreement between the modified ST and the iV-body results is remarkable. On all 
measurable scales it is better than 2-3%. This is important for analysing future cluster 
surveys because it means that existing semi-analytic methods can be used instead of 
having to perform time consuming simulations for all neutrino masses. 

Alternatively, neglecting the neutrino perturbations in the A^-body simulation will 
also be a very accurate approximation for '^m,, < 0.5 eV as long as only matter halo 
properties are considered. This approximation is not valid for a precise calculation 
of the matter power spectrum [10]. In general the accuracy of the approximation is 
determined by contrasting the neutrino free-streaming length with the physical extent 
of the scales simulated: Considering halo properties and realistic neutrino masses, this 
approximation is very good. 
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